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Abstract 



Accurate (^$ 1%) predictions for the anisotropy of the Cosmic Background 
Radiation (CBR) are essential for using future high-resolution (< 1°) CBR 
maps to test cosmological models. In many inflationary models the variation 
("running") of the spectral index of the spectrum of density perturbations is a 
significant effect and leads to changes of around 1% to 10% in the CBR power 
spectrum. We propose a general method for taking running into account 
which uses the derivative of the spectral index {dn/d In k). Conversely, high- 
resolution CBR maps may be able to determine dn/dlnk, giving important 
information about the inflationary potential. 
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Typeset using REVTeX 



The cosmic background radiation contains a wealth of information about the spec- 
trum of primeval density perturbations. This is because CBR anisotropy on a given 
angular scale arises largely due to density perturbations on a (comoving) length scale 
L ~ (0/l°)lOO/i _1 Mpc. Since the COBE detection of CBR anisotropy on angular scales 
of 10° to 90° |IJ, more than ten additional detections on angular scales from about 0.5° to 
20° have been reported @. In addition, plans are being made for a satellite-borne exper- 
iment within the decade that will map the CBR sky with an angular resolution of better 
than 1° and an accuracy that is an order of magnitude better than current measurements 
H . Thus, in the near future CBR anisotropy should be able probe inhomogeneity on length 
scales from about Mpc to 30, OOO/^ 1 Mpc. 

A key to using CBR measurements to reveal the underlying spectrum of density per- 
turbations is the accurate calculation of the expected anisotropy in a given model. Much 
progress has been made in understanding and taking into account all the relevant micro- 
physics ||, and several groups are now making a concerted effort to calculate expected CBR 
anisotropies with an accuracy of better than 1% ||. 

Much of this effort is directed at inflation, as CBR anisotropy has the potential to 
both test the inflation hypothesis and reveal important information about the underlying 
scalar-field potential Inflationary models predict approximately scale- invariant spectra 
of density (scalar metric) perturbations [|]] and gravity-wave (tensor metric) perturbations 
||, and both contribute to CBR anisotropy. The following parameters have been identified 
as important for accurately computing the expected anisotropy ||: the power-law indices 
of the scalar and tensor spectra, tit ~ and n « 1; the overall amplitudes of the scalar 
and tensor perturbations, often quantified by their contributions to the variance of the 
quadrupole anisotropy, Qs and Qt] the Hubble parameter h = if /100kms _1 Mpc -1 ; the 
baryon density, which is constrained by primordial nucleosynthesis to the interval Qsh 2 — 
0.009 — 0.022 [|Kj; and possible contribution of a cosmological constant to the energy density 
of the Universe today fl^- (In addition, some have considered the possibility of a total energy 
density less than the critical density predicted by almost all models of inflation, nonstandard 
ionization histories for the Universe, and variations in the nonbaryonic component of the 
matter density, e.g., adding a small amount of hot dark matter.) 

In this paper we emphasize that the spectral indices n and in general vary with 
scale and point out that for many interesting models of inflation (chaotic, natural, and 
new) the variation in scalar spectral index leads to significant corrections (1% to 10%) in 
the predicted CBR anisotropy. Conversely, this means that a high-resolution CBR map 
could be used to extract information about the variation of n with scale and thereby reveal 
additional information about the inflationary potential. We thus make the case that the 
variation of the scalar spectral index should be taken into account when calculating CBR 
anisotropy, and suggest that it is most sensibly done by using dn/d Ink. 

CBR anisotropy on the sky is usually expanded in spherical harmonics, 

5T(n)/T = J2ai m Y lm (n). (i) 

Inflation makes predictions about the statistical properties of the multipole moments; since 
isotropy in the mean guarantees that (a/ m ) = and the underlying perturbations in almost 
all inflationary models are gaussian, the variance C[ = (|a/ m | 2 ) serves to specify all statis- 
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tical properties. (Here and throughout brackets refer to the average over an ensemble of 
observers.) Measurements of the CBR temperature on the sky can be used to estimate the 
statistical properties of the underlying density perturbations. In particular, the Cj's can be 
estimated. Because the sky is but a finite sample, a fundamental limit to the accuracy of 
the estimate (referred to as cosmic variance) is given by 



or 2 

u i)i 2l + 1 



(2) 



Other major (and presently dominant) sources of uncertainty include receiver noise, various 
instrumental systematic errors, foreground sources (our own galaxy, radio sources, etc.), 
limited sky coverage, and finite resolution (a map with angular resolution 9 is only sensitive 
to multipoles with / < 200°/#). 

High-resolution maps of the CBR probably offer the best means of studying the scalar 



and tensor metric perturbations predicted by inflation [pTT| , p~2[| . Such maps may also provide 
valuable information about the Hubble constant Hq, the cosmological constant A, the baryon 
density Qb, and the total density of the Universe Q; however, other measurements will 
complement the determination of these parameters. If the four parameters describing the 
scalar and tensor perturbations are measured to some level of accuracy, properties of the 



underlying inflationary potential V((p) can be determined [O] 



Vn = 1.65Qr m-p\ 



Vx = ±V-8Tcn T V N /m P i = ±J V N /m Ph 



= 4ir[(n - 1) - 3n T ] V N /m Pl 2 = 4tt 



[n - 1) + -r 



Vjv/rapi 2 



(3) 
(4) 
(5) 



where r = Qt/Qs-, a prime indicates derivative with respect to <p, and the sign of V is 
indeterminate. In addition, a consistency relation = —r/7 must be satisfied, and the 
factors of | arise from using it |14[] . Subscript N indicates that the potential is to be evaluated 
at the value of where the scale corresponding to the present Hubble scale {k^ = Hq) crossed 
outside the horizon during inflation. This generally occurs around iV ~ 50 e-foldings before 
the end of inflation, though the precise expression depends upon the model of inflation, the 
reheat temperature, and any entropy production after inflation. (Only the expression for 
depends upon the definition of iV; the other two always apply.) The expression for iV 
can be written as 



N ~ 54 + - ln(-n r ) + - ln(T RH /10 6 GeV) - - In 7 - In h, 
6 3 3 



(6) 



where Trh is the reheat temperature, 7 is the ratio of the entropy per comoving volume 
today to that after reheating which quantifies any post-inflation entropy production, and 
the perturbation spectrum has been normalized to COBE. (In calculating iV it has been 
assumed that inflation is followed immediately by a matter-dominated epoch associated 
with coherent oscillations of the inflaton field and then by reheating.) 

The above expressions were derived in a systematic approximation scheme that relates 
the derivatives of an arbitrary smooth inflationary potential to CBR observables [i"3"|,|H) . 
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The expansion parameter is the deviation from scale invariance, and formally involves all 
the derivatives of the potential, m-p\ n V^ /V^ (constant V corresponds to the scale-invariant 



limit; see Ref. |L6[] for a discussion of this scheme). For most potentials the deviation from 
scale invariance of the scalar and tensor spectra, quantified by (n — 1) and n?, serve as the 
expansion parameters. The above expressions are given to lowest order in (n — 1) and Ut', 
the next-order corrections are given in Ref. 
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The crucial point for the present discussion is that the spectra of scalar and tensor per- 
turbations are only exactly power laws for an exponential potential. In general, they vary 
with scale, though dn/dlnk and drix/dlnk are second order in the deviation from scale 
invariance, i.e., involve terms that are 0[(n — l) 2 ,n^, (n — l)nj-]. Since the present data 
indicate that scalar perturbations do not differ from scale invariance by a large amount, 
n — 1 = 0.10 ± 0.32 |T7j, the variation of the spectral indices is expected to be small. Fur- 
ther, indications are that the tensor perturbations are subdominant and in any case only 
contribute significantly to multipoles I — 2 to 50. However, we shall show that the vari- 
ation of scalar spectral index is important, given the desired precision for the theoretical 
predictions of the multipoles. 

The power spectrum for the scalar perturbations is given by 



P(k) = Ak n{ - k) 




n.jv + ln(k/k]y)(dn/d In k) - 



(7) 



The contribution to the Ith multipole comes from wavenumbers k centered around Z/to, 
where r ~ 2 /H Q is the distance to the last scattering surface. Recalling that the char- 
acteristic scale fcjv w as chosen to correspond to the current horizon size, this implies an 
approximate scaling relation for the C/'s which relates them to a spectrum with constant 
spectral index: 

/ , \ \n(l/2)dn/d\nk 

CMk)]"!-) Q[n(k)=n N }. (8) 



If \dn/d\nk\ >3 x 10~ 4 , the effect of ignoring the "running" of the spectral index over the 
range 1 — 2— 1000 is greater than one percent, which is significant compared to the accuracy 
goal for CBR anisotropy 0. We now show that values this large are expected in interesting 
inflationary models. 

In general, the derivatives of the scalar and tensor spectral indices are related to the 
inflationary potential and its derivatives. The lowest-order expression for all the derivatives 
of n and tit can be obtained by simply differentiating the lowest-order expressions, 



<t>=<t>N 



n — 1 = nx + 

using the fact that to lowest order 

d 

dliak 



mpi d ( mp\V r 



V v 



(9) 



(10) 



4>=<j) N 




4>=4>n 
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If n and tit are expressed as a function of N, one can use the fact that d/dlnk = —d/dN 
to obtain the desired derivatives even more easily. 

It is thus a simple matter to obtain the first derivatives of n and n T : 



dn 



m Pl 3 V"> ^ f mpiV' ^ _L_ ( m Pl 2 V" ^ / m P1 \/' x 2 



dink 32tt 2 V V J \ V J 8n 2 \ V J \ V 

driT 



1 f m Pl 2 V" \ ( m PX V' \ 1 ( m PV V' \ 
V J [ V J ~ 32^ { V J 




dink 32tt 2 

Equivalent expressions can be obtained by using the previous equations relating r and n — 1 
to the potential and its first two derivatives: 

dn 1 /27T /m P1 V;\ /- 4. 6 2 

*i= T iw Thrr + 7 ( ""- 1)r+ s r ' (13) 



AT 

rfluT = - nr[(n - 1) - nr] = 7 



x 1 ■ 

(n - 1) - -r 



(14) 



where the upper sign applies if > and the lower if < 0, and the factors of | arise 
from using the consistency relation rir = —r/7. From these expressions we see that the size 
of both dn/dliak and drix/dlnk is controlled by the ratio of tensor to scalar perturbations, 
and further, that the size of drix/dlnk depends upon the difference between n — 1 and n*r, 
which in many models is small. 

We now quantify expectations in several popular models of inflation. As noted earlier, 
for an exponential potential dn/dlnk = dn^jdXYyk = 0. For inflation models that are based 
upon Coleman- Weinberg like potentials, V{(j>) = Ba 4 /2 + J B0 4 [ln(0 2 /a 2 ) -1/2], 

dn/dlnk ~ -1.2 x 10" 3 (50/iV) 2 , (15) 
d m n/d\nk m ~ -3m!/A^ m+1 . (16) 

Chaotic-inflation models are usually based upon potentials of the form V((f>) = a<p b (a is a 
constant and b = 2, 4, • • • is an even integer) and 

dn/dlnk = -4 x 1(T 4 (6/2 + l)(50/iV) 2 , (17) 
d m n/d In A;" 1 = -m! (6/2 + l)/iV m+1 . (18) 

For the interesting cases of b = 2 and 4, dn/dlnk = —0.8 x 10~ 3 (b = 2, N = 50) and 
— 1.2 x 10~ 3 (6 = 4, iV = 50). Finally, for the "natural" inflation model, where V(4>) = 
A 4 [l + cos(0//)], the following approximate expression applies for / < mpi (which is the 
regime where the deviation from scale invariance is significant and 1 — n ~ mpi 2 /8nf 2 ): 

7T 2 

rfn/rflnA; = — ^-(« ~ l) 2 exp[iV(n - 1)]. (19) 

Varying (1 — n) from 0.04 to 0.3 and N from 40 to 50, dn/dlnk varies from about — 10~ 7 
to almost — 10~ 3 . Even though 1 — n can be large in these models, r is very small when it 
is, and dn/dlnk never approaches (1 — n) 2 . 
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It is not a complete surprise that dn/d In k is similar in all these models. On naive grounds 
one might expect that (n — 1) oc 1/N m , so that dn/dlnk = —dn/dN = m(n — 1)/N. This 
is true for new and chaotic inflation where m = 1. (As usual, the situation with "natural" 
inflation is more complicated.) 

Are there models where running is more important? For an ad hoc potential the answer 
is yes; consider V(<f>) = Vbexp(— acf) b ). Here 



(n-l) -- 



(2 - by 



ab 
8^ 



2-a 



1 a 



where a = 2(1 — 6)/(2 — b). For 6 ^ 1, n — 1 can be large and dn/dhxk ~ a(n — 1)/N. 

As mentioned earlier, the running of the tensor spectral index is expected to be less 
important because the tensor perturbations are likely to be subdominant and only contribute 
significantly for / < 50; in addition, dnx/dlnk is smaller (being proportional to the difference 
between n — 1 and tit which is often small). For the potentials discussed above dnx/dlnk = 
(exponential), -2 x 10- 7 ((x/ m Pi) 4 (50/iV) 4 (new), -2.0 x 10~ 4 6(50/A^) 2 (chaotic), and 
-7r 2 /4(n - l) 2 exp[iV(n - 1)] ("natural"). 

Figure 1 displays the CBR angular power spectrum for 6 = 6 chaotic inflation (where 
dn/dlnk = —0.0020), calculated without and with the running of the scalar spectral index. 
The correction due to the running of scalar index is significant (about 10%) and potentially 
measurable. The results shown have been calculated using the power spectrum in Eq. (^). 
We have also calculated the C/'s using the approximation in Eq. (§), and the maximum error 
in any C\ is less than 0.6%. Thus, for applications requiring 0(1%) accuracy, it should be 
sufficient to calculate a model with a fixed n and then scale the results according to Eq. ([8]) 
to obtain C/'s for dn/dlnk ^ 0. We also note that the correct fc-space power spectrum is 
simple to include in any Boltzmann code. 

In summary, expectations for \dn/dlnk\ in popular inflationary models range from —2 x 
10~ 3 to around —4 x 10 -4 . Of course, the value of dn/dlnk in "the model of inflation" could 
be larger or smaller. At the high end of this range, neglecting the running of scalar spectral 
index leads to errors of 10%, more than an order of magnitude larger than the accuracy 
desired The running of the scalar spectral index can be into account easily, accurately, 
and with generality by using dn/dlnk. Based upon the models we have looked at one could 
adopt dn/dlnk ~ (n — 1)/N as a default estimate. 

If the running of the scalar spectral index is large enough to detect, the third derivative 
of the scalar potential can be measured |L6|| : 



Vjf/mpi = ±39v^ [-7{dn/d In k)/r + 0.9r + 4(n - 1)] Q T . (20) 

The feasibility of determining dn/dlnk from a high-resolution map of the CBR sky is cur- 
rently under study [15]. 

Two final points. First, what about the next-order corrections? They involve 0[(n — 
l) 3 , riy, • ■ ■] terms: corrections to dn/dlnk, n, Qt, Qs and the d 2 n/dlnk 2 term in the 
expansion for n. Provided that the deviation from scale invariance is not too large, they 
should be small (less than about 1%) because they are suppressed by an additional factor 
of 0[(n — l),nx]; e.g., dn 2 /dink 2 = —5 x 10~ 5 (50/iV) 3 for new inflation, which leads to a 
correction at / = 1000 of about 0.5%. If the d 2 n/dlnk 2 should be larger, its size might be 
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turned to good purpose; because of the qualitative difference between it and the d/d\nk 
term it might possibly be measured, revealing information about the fourth derivative of 

Last, but perhaps not least, the running of the scalar spectral index is also of some 
relevance when extrapolating a COBE-normalized spectrum to astrophysical scales; e.g., 
the correction to as is about —3% for dn/dhxk = — 10~ 3 [EU|. 
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FIG. 1. Predicted angular power spectra for 6 = 6 chaotic inflation with (solid) and with- 
out (broken) the running of the scalar spectral index (n = 0.92, h = 0.7, Qb = 0.025, and 
dn/dlnk = -0.002). 
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